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ABSTRACT 



We have calculated the opacities and secondary production mechanisms of 
high energy photons arising in gamma-ray burst internal shocks, using exact 
cross-sections for the relevant processes. We find that for reasonable choices 
of parameters, photons in the range of lO's to lOO's of GeV may be emitted in 
the prompt phase. Photons above this range are subject to electron-positron pair 
production with fireball photons and would be absent from the spectrum escaping 
the gamma-ray burst. We find that, in such cases, the fireball becomes optically 
thin again at ultra-high energies (> PeV). On the other hand, for sufficiently 
large fireball bulk Lorentz factors, the fireball is optically thin at all energies. 
Both for 77 self-absorbed and optically thin cases, the escaping high energy pho- 
tons can interact with infra-red and microwave background photons to produce 
delayed secondary photons in the GeV-TeV range. These may be observable with 
GLAST, or at low redshifts with ground-based air Cherenkov telescopes. Detec- 
tion of the primary prompt spectrum constrains the bulk Lorentz factor, while 
detection of delayed secondary gamma-rays would provide a consistency check 
for the primary spectrum and the bulk Lorentz factor as well as constraints on 
the intergalactic magnetic field strength. 

Subject headings: gamma rays: bursts — gamma rays: theory — radiation mecha- 
nisms: non thermal 



1. Introduction 



GeV photons have been detected from a number of gamma-ray bursts (GRB) (Dingus 
2003) with the EGRET detector on the Compton Gamma Ray Observatory. Photons at 
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energies in this range and higher can make pairs via 77 — > e+e~ interactions in the GRB 
fireball, the optical depth being dependent on the bulk Lorentz factor r^, so spectral obser- 
vations at this energy can provide useful constraints (Baring & Harding 1997; Lithwick & 
Sari 2001) on this key ingredient of GRB models. In addition, high energy photons which 
escape 77 absorption in the source can still suffer such interactions far from the source 
but before reaching the observer, against cosmic infrared background (GIB) and/or cos- 
mic microwave background (GMB) photons. The electron-positron pairs thus produced can 
inverse- Gompton scatter again on the most numerous GMB photons, giving rise to a delayed 
secondary spectrum (Plaga 1995; Dai & Lu 2002), which ranges typically up to lOO's of GeV. 

In this paper we investigate, using the exact cross sections, the high energy interactions 
of photons produced in a simple GRB internal shock model, as a function of the luminosity, 
the typical variability timescale and the bulk Lorentz factor. We approximate the original 
GRB spectrum through an observationally motivated Band function (Band et al. 1993), 
and calculate the emerging spectrum as a function of the burst parameters. This spectrum 
is then subjected to interactions with the GIB and GMB, and the shape of the secondary 
spectrum and its time delay is calculated assuming typical cosmological distances ^ ~ 1 and 
various values for the intergalactic magnetic field strength. The pair-production cut-off of 
the primary spectrum and the delayed secondary spectrum fall, for typical GRB parameters, 
in the GLAST (Gehrels & Michelson 1999) energy range, and for a range of intergalactic 
magnetic field strengths the latter should be detectable by GLAST as well as by the newer 
generation of air-Gherenkov telescopes (Weekes 2001). 



2. GRB internal shock model 

Denoting with U the GRB fireball total energy density, after the initial expansion phase 
the kinetic energy is carried by the baryons since rup ^ irie, where nip and nie are baryon 
and lepton masses respectively. Thus, U Up An'pnipC^Vl, where n'^ is the total baryon 
number density in the comoving fireball frame and is the bulk Lorentz factor of the fireball 
in the observer's frame. ^ Later on, a significant fraction of the fireball's kinetic energy is 
transferred to leptons, in internal shocks which randomize the relative bulk motion between 
different portions of the expanding fireball, as well as external shocks when the fireball is 
decelerated by the external medium. Here we concentrate on the internal shocks, which 
are thought to be responsible for the usual prompt MeV emission. The fraction of energy 



^We use primed variables in the comoving frame and un-primed variables in the laborotary or observer's 
frame 
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transferred to the leptons is given by a parameter — Ue/U ^ rrie/mp. Charged particles, 
p and e"'"e~, are accelerated in the shocks by interaction with magnetic fields within the GRB 
fireball. High energy protons and electrons cool mostly due to synchrotron radiation, but the 
proton cooling time is much longer than that of the electrons, and the prompt 7-ray emission 
is most likely due to electron synchrotron emission, whose cooling time in the prompt phase 
(tens of seconds after the trigger) is generally shorter than the dynamical or expansion time. 

In this so-called fast cooling case, valid in the internal shocks, most of the electron energy 
is converted into 7-rays: ^ = L^^iso/ (4vrrgj^c). Here = 2c6trl is the shock radius, 6t 
is the observed variability time scale, and L^,,iso is the isotropic-equivalent 7-ray luminosity in 
the observer's frame during the prompt 7-ray emission phase. Neglecting pair formation 
in the shocks, the total volume number density of leptons and baryons in the fireball (which 
are the same due to charge neutrality) can be calculated, assuming Up ~ CZ-y/Sg, as 

in the comoving frame. Here 7p ~ 1 is the comoving random Lorentz factor acquired by 
protons in the internal shocks. The magnetic field energy density Ub = -B^/Stt can be 
estimated as 



so in the comoving frame Eb — Ub/U is the fraction of the shocked fireball energy that goes 
into magnetic energy. 

The process of Fermi acceleration in the shocks leads to a differential electron and proton 
volume number density distribution which is a power law in energy, 

ie,p /e,p;min \ /e,p;min/ 

where 



^ J £e7;,min ; T^.min ~ 1 (4) 



are the minimum Lorentz factors for electrons and protons respectively in the comoving 
frame. The maximum Lorentz factors are obtained by equating the acceleration time {t'^ ^ — 
Arrijjj / eB') of a particle j to the shorter of the synchrotron cooling time 

^ 67rm^ 
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and the dynamic time {StFj,) in the comoving frame. Electrons cool faster, by synchrotron 
radiation, than the dynamic time and the maximum electron Lorentz factor is 



f I 67re 



7e,max y j^^^^^f 

^ 2.2 X 10^^-/^<//.;/^,e-Zt5ti/^r,^/i (6) 



in the comoving frame. Here we have chosen L^^iso = 10^^ ^^^52 erg/s, = 10 ^£e-i, 
Eb = IQ-'^Eb-i, St = l0-^5t^2 s, A = lOAi and Tb = lO^-^Ffe^s.s. Note that the inverse 
Compton (IC) cooling time of the electrons is proportional to the synchrotron cooling time 
in equation (5) and the proprtionality constant, in the fast cooling case, is given by F = 
(— 1+ 1/1 + 4£e/£B)/2 (Sari & Esin 2001). For the present choice of parameters, Y ki 0.6 and 
we have ignored the IC coohng of the electrons compared to the synchrotron coohng while 
calculating 7e,inax- Our choice of the value of £5 is motivated by the assumption that the 
observed peak 7-ray energy is due to electron synchrotron radiation and only a small fraction 
of the total 7-ray power is radiated by IC process at small internal shock radii (Daigne & 
Mochkovitch 1998; Derishev, Kocharovsky & Kocharovsky 2001; Baring & Braby 2004). A 
recent study by Zhang, Kobayashi & Mcszaros (2003) also suggest that the fireball value of 
Sb should be 10-100 times larger than the afterglow fit of 0.01-0.001. 

The observed 7-ray spectrum from the GRB prompt emission can be approximated by 
the phenomenological Band function in the high energy range as ~ e~° above the synchrotron 
peak energy e^^pk, which is typically a few hundred keV. For observed cases of a > 2, most 
of the 7-ray energy is concentrated near e^,pk and one can roughly estimate the peak volume 
number density of photons in the fireball as 

^ 47rr2i^crbe^,pk 

in the comoving frame. Below the synchrotron peak energy the Band spectrum is ~ and 
the total number density of photons in equation (7) increases only logarithmically. 

Low energy photons are expected to be absorbed by fireball electrons in the presence of 
a magnetic field by the synchrotron self- absorption mechanism below a photon energy of 
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in the observer's frame. Here we have estimated e-y,ssa for power law indices in equation (3) 
of p = 2,2.5,3 (top to bottom in the equation above). Observed GRB spectra, fitted by a 
Band function, do not correspond in a straightforward manner to a theoretical synchrotron 
spectrum from a one-zone region with some power-law index p, especially in the fast-cooling 
regime. However, the usual Band function values can arise as a result of the superposition of 
many shocks, each giving rise to a power law energy distribution with index 2 < p < 3. In our 
calculation, we assume the nominal Band fit values for the resulting photon spectrum, and we 
use a nominal reference value for the absorption energy of e^^ssa — lOe^^i keV, which is more 
than an order of magnitude below the observed synchrotron peak energy of e-y,pk = 500e-y,2.7 
keV. 



3. High energy photon interactions and internal attenuation 

Observation of high energy photons from the GRB fireball is mainly limited by the 
opacity of two photon annihilation into an electron-positron pair (e^). The cross-section for 
77 — > in the non-relativistic (NR) and extremely relativistic (ER) cases is given by Jauch 
& Rohrlich (1955) as 



a';^(u;) ^ -(7Thi ' ' 



CH = J^Th^ ln£^-l) (9) 



mlc^ / 2uj 

„- ill 5~ I ill n 

8 o;^ \ meC^ 

where u is the photon energy in the center of mass frame. The observed photon energy 
{E.y) above the threshold for pair production with photons at an energy must satisfy the 
condition uj > 171^.0^ or > -E^^th = ''^eC^r^/2e^. Observation of photons of a maximum 
energy e^^max > -£'7,th may therefore be used to put a limit on if the corresponding 
optical depth is unity. Photons below E^^^^, may escape from the fireball (provided no other 
interaction becomes important) even though they might interact with photons from the high 
energy tail of the Band function, but the relevant optical depth is less than unity because of 
a much lower photon density, at that energy, than in equation (7). 

For a given F^ and 5t, the lower limit threshold of the 77 absorption energy range within 
the source can be found from the pair production threshold energy with peak synchrotron 
photons of energy e^^pk as 

m^c^F^ 

£;7,pk,th = ^ 26F^,2.5S,i7 GeV (10) 

in the observer's frame, for the present choice of parameters. The total volume number 
density of photons and hence the 77 optical depth in the fireball increases slowly as the 
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photon energy decreases from e-y^pk to e-y,ssa- High energy photons in the energy range capable 
of producing pairs with fireball photons in the energy range e^^pk - e7,ssa ^-re absent in the 
spectrum escaping from the source. Ultra high energy photons, however, may escape the 
GRB fireball as the pair production cross-section, in the ER limit in equation (9), decreases 
with increasing photon energy and so does the optical depth. Ultra high energy primary 
photons will appear again above this "thinning" energy. This energy can be found, roughly, 
from the optical depth corresponding to the 77 cross-section with photons at self absorption 
energy (e-Y,ssa) in the ER limit in equation (9) as 

_ 3A L^^isoQ-Thmlc^ 

« 2 X 10%L^,52r,-i55t:^e7j GeV (11) 

where A = ln(2 2e^^ssaE-y/ meTb) — 1 ~ 3A3, for the same parameters used in equation (11). 

We calculate the high energy photon interactions with other fireball photons and elec- 
trons using the exact cross-section formulae (Jauch & Rohrhch 1955). The opacities of high 
energy photons are plotted in Figure 1, as a function of observed energy {E^), due to electron 
Compton scattering, pair production with electrons (67 — > e^) and by 77 interactions, 
dominantly, with fireball photons in the energy range e^^ssa - The cross-section for the 
67 — > process is given by Jauch & Rohrlich (1955) as 



<7,7(£;;) = 2.25 X 10-3^<7Th ( 



in the non-relativistic and extreme relativistic limits respectively. Here E'^ is the photon 
energy in the comoving frame. The fiat region, in Figure 1, for the 77 pair creation curves 
above £'-y,pk,th corresponds to a change in the target photon energy from e-y,pk to e-y,ssa and a 
slow increase of the total target photon number density n'^, in equation (7), accordingly. 

We have plotted the boundaries of the optically thin and thick regions for high energy 
photon emission from GRB internal shocks in Figure 2, for different choice of parameters 
L-y-i e-y,pk and 5t. Note that for certain parameter values, e.g. L^^i^o — 10^^ erg/s, 5t — 0.1 s 
and e^^pk = 0.5 MeV; or -L^,iso = 10^^ erg/s, 5t — 1 s and e-y,pk = 1 MeV, the GRB fireball is 
optically thin to photons of all energies above 775 or 676 respectively. 

Extra complications can arise due to increased pairs produced in the fireball by 77 
interactions by high energy photons (Totani 1999; Meszaros, Ramirez- Ruiz & Rees 2001), 
which can affect the photon opacities plotted in Figure 1. This is a complex problem. 
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requiring a numerical treatment (Pe'er & Waxman 2003). Here we use an approximate 
treatment, using a pair formation rate in the comoving frame of ca^j{e'^)dN!^ J de'^ to estimate 
the total number of pairs created by an incident (denoted by "i") high energy photon. 
The differential volume number density of incident high energy photons capable of producing 
pairs, as a function of energy, can be estimated using the Band spectrum and equation (7) 
as 

^ ~ < - e' > e' (13) 

7,pk IS/ 7,pk^ 

Here e!y,pk,th comoving threshold energy for pair production which is similar to the 

threshold energy in equation (10) in the case of synchrotron peak photons as targets. The 

target (denoted by "t") photon number density, in general, may have a different energy 
distribution dN'^^^/dE'^. The total volume number density of pairs created in the fireball is 
thus 



where e^^, j^g^ ~ 7e,max"^eC^ is defined in equation (6), and e'^^^-^ = mlc^/2E'^. For this estimate, 
we approximate the 77 cross-section as a^^ ~ (3/8)(TTh at pair production threshold 
with peak synchrotron photons (e-y,pk). This leads to a simplified expression for the total 
volume number density of pairs created in the fireball by integrating over the photon energy 



distribution above e'^^^^, 



The total number density of target photons in equation (7) increases by a factor ~ ln(e^^pk/ e7,ssa) 
as the target photon energy decreases from e-y^pk to e-y,ssa in case of a = 2. However, one 
needs incident photons of energy much higher than e^^pk to produce pairs with photons 
at e^.ssaj and the incident number decreases oc Thus equation (15) is estimated for 

targets which are at threshold with incident photons of energy > e^^pk- We can express the 
secondary pair density as the ratio to the initial electron density in the fireball, 

/or / 9^2 \ a-l 

^ 6 (7Th -^7,iso^e^p / ^S.pk 



1287ra-l TlStel^^ 

786L^,52£e,-i<^t^r,j5; a = 2 

2.3L7,52£e,-i^t:2e7,2.7r5j.5; Qi = 2.5 (16) 
0.008L7,52£e,-i5C^e2 2.7r-2'.5; a = 3 
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using equation (1). 

Shock accelerated high energy electrons may also produce 6=*= pairs interacting with 
synchrotron photons (e7 ee^). (The cross-section for this process is given by equation 
(12) with the replacement E'^/nie — > ^'^E'^/me). However, the number of secondary pairs 
produced in this case is negligible compared to the initial lepton number density because of 
the lower number density of high energy incident electrons. 

The secondary pairs can annihilate with themselves or with the leptons originally carried 
in the fireball, e"^e~ 77, if the number density of pairs increases substantially. The cross- 
section for this process is (Jauch & Rohrlich 1955) 



8 7^(7: + 1) 



(7^ + 4 + 1) In (7: + ^^?^) - (7; + 3) 



(17) 



Because of a decreasing cross-section with energy, low energy pairs annihilate faster. 
The annihilation time for an incident electron can be estimated as t'^^^ = l/(ccre±n'_j_) in 
the comoving frame. However, synchrotron cooling is even faster and the ratio of pair 
annihilation time and synchrotron cooling time, for example at 7^ = 5 (where the electron is 
still relativistic and synchrotron loss formula applies), Ff, = 10^'^, St = 0.01 s and in a = 2 
case in equation (16), is te,syn/^ann ~ 0-1 using equations (5) and (17). This ratio decreases 
at higher energy. Hence electrons lose most of their energy to synchrotron photons within 
a dynamical time StTh in the comoving frame. As a result, secondary pairs produced in the 
fireball cool down to nonrelativistic values and drift along roughly with the bulk Lorentz 
factor r^. We consider here situations where the secondary pair density, equation (16), is at 
most comparable to the original fireball lepton density, in which case the scattering depth 
is not substantially changed, and any photons from annihilation of cold pairs lead to higher 
generation pairs whose number similarly does not affect the scattering optical depth. We 
note that in the Thomson limit, inverse Compton (IC) cooling time is of the same order as 
the synchrotron cooling time t'ejc/^e.syn ~ ^e/^B with an increasing value at higher energy 
and in the Klein- Nishina limit. 



4. High energy photon spectrum and interaction with the diffuse background 

Using the phenomenological Band function spectrum for the intrinsic high energy photon 
flux from the GRB, the spectrum that would reach the observer after taking into account 
in-source 77 attenuation at a luminosity distance Dl is 

(PN^ {a - l)L-^,iso / e 



(^) e--— (^-); 6,>6^,pk (18) 
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where Tyy^grhi^y) is P^i^ production optical depth (see in Figure 1) in the GRB fireball. 

The high optical depth in the fireball, however, reduces substantially the ultra-high energy 
photon flux above the energy -E^,pk,th given in equation (10) and below the energy -£/-y^thin 
given in equation (11), while for some parameters the GRB is in the region of Figure 2 
where the escaping spectrum is optically thin at all energies. Here we consider high energy 
photon emission from the GRB fireball for two cases. One is an example where the escaping 
spectrum is optically thin at all energies (e.g. L^ igo = 10^^ erg/s, St = 1 s, e^^pk = 1 MeV 
and Fb = 10^-^ in Figure 2). The other example is a case when the fireball is opaque to 
photons in an energy range between £'7,pk,th and £'-y,thin (e.g. iv-y,iso = 10^^ erg/s, 5t — 0.1 s 
and e^,pk = 0.5 MeV and F^ = lO^-^). 

For typical GRB locations at redshift z — 1, most high energy photons above ~ 70 GeV 
produce pairs (Salamon & Stecker 1998), reaching an energy up to ~ ?7T.e7e,max/2, due 
to 77 interactions with cosmic infra-red background (GIB) photons (in the ~ 100 GeV-TeV 
range), and cosmic microwave background (GMB) photons (in the PcV range). These high 
energy pairs, in turn, cause delayed high energy photon emission by inverse Compton (IC) 
scattering of GMB photons (Dai & Lu 2002). Here we follow the treatment of Dai et al. 
(2002) to calculate the delayed high energy photon spectrum in the GeV-TeV range. 

The corresponding time integrated flux of electrons (and positrons) from high energy 
GRB photons interacting with GIB and GMB photons is 



d7e 27: Dl (2me)"-i 



X [1 - e-^-^-^(^>-)] e ( < 7e < ^) , (19) 

where is a step function. We have assumed that each electron and positron of the pair 
share 1/2 the photon energy, i.e. 7e = e-y/2me. We calculate the pair production optical 
depth with background photons (GIB and GMB) at threshold (2e^e^^cib;cmb — ^1) as 

T77,bkg(e7) = max{ncib(mg/2e^,cib),?T'cmb("^e/2e7,cmb)} x (3/8)(7Th-DL (20) 

The total number density of GIB photons is fitted from the intensity 1 (W-m~^ sr~^) given 
in Malkan & Stecker (2001) as 

47rJ(l + zf 

?T'cib(e7,cib) = ■ (21) 

C£7,cib 

The GMB photons, on the other hand, have a black body distribution peaking at an energy 
of 2.7(1 -|- ^) K and a total volume number density oc {1 + z)^. 
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The secondary electrons in equation (19) are cooled by IC scattering on CMB photons 

on a timescale tic(7e) = 3mec/(47e(TTh'Wcmb) ~ 7.3 x 10^^(7^/10^)"^ s in the local rest frame 
(Dai et al. 2002), where ticmb is the local black body CMB photon energy density. The 
total integrated, first generation, IC photon spectrum in the observer's frame, convolving 
the electron spectrum with the IC spectrum from a single electron in the Thomson limit 
(Blumenthal & Gould 1970), is then 



rf'A^delayed IC _ f f ^^dN, £ 



dtdE^ J J d'ye dtdE^ 

3(7ThC (« - l)<pk 



32nDl (2me)"-^ 



X 



-^7,iso^grb CXp [ T^^^bkg (-E^)] 



J d-fe% ^""^^^ exp [-T^^,grb(27eme)] (1 - exp [-r^7,bkg(2me7e)]) 



At(7, 

^ J ^cmb(e7,cmb) [r, ZT- 1„ "^7 , r i /l-,2, ^7 



^cmb l,^7,cmb ) 



7,cmb 2 



7,cmb 



2£;^ln-^^^ + E^ + A-il 

^7,cmb 



C7,cmb 

,e ..... 27 

e C7,cmb 



2) 



Here the integration over 7e ranges from 7e,min = max{me/(4e-y,cmb), s/ E.y/ e^^caih/'^} to ~ 
7e,max/2 whcrc 7e,max IS the maximum electron energy in equation (6). The maximum time 
scale for delayed emission is At(7e) = maxjAtjc, At^, At^, tgrb}, where Atic(7e) = 'fic/(27e) 
is the IC cooling time; At^(7e) = A^^^bkg/(27pC) is the angular spreading time; AtB{je) = 
^ic^b/2 is the delay time due to magnetic deflection, all evaluated in the observer's frame; 
and A^^^bkg = (niax{ncib, 'n'cmb}(3/8)(TTh)^^ is the pair production mean free path against 
background photons. For 10~^°i?iG,-2o G inter galactic (IG) magnetic field, the defiection 
angle, for 9b < 1, is 9b ^ K/Rg ~ 1-3 x 10-^(7e/10^)-2(5iG -2o)^ where = -femeC^/ieB) 
is the Larmor gyroradius of the electron. We have plotted different time scales as a function 
of secondary pair Lorentz factor (7e) in Figure 3 for tgrb = 50 s and several IG magnetic 
field values, at = 10^^ cm corresponding to redshift 2; ~ 1 in an fl^ — 0.3, fl\ — 0.7 and 
Ho = 75 km s~^ Mpc~^ cosmology. 

We used two parameter sets {Lj^iso — 10^^ erg/s, 6t — 1 s, e-y,pk = 1 MeV and = 
10^-9) and (L^,iso = 10^^ gj^g/g^ ^ g.l s and e^,pk = 0.5 MeV and = lO^-^). These 
correspond to a case when the primary spectrum escaping from the GRB fireball is optically 
thin to photons of all energes, and a case when the primary fireball escaping spectrum 
is optically thick to 77 against its own photons in a band of photon energies, shown in 
Figures 4 and 5 respectively. The prompt primary spectra, using equation (18) multiplied 
by exp[— r^^^bkg(£7)] to allow for attenuation in the background radiation fields, are plotted 
as the dark sohd curves in both figures. We have also plotted the first generation delayed 
spectra from IC scattering of secondary leptons on the CMB, by numerically integrating 
equation (22) for two different IGM magnetic field strengths (Bjg) of 10"^° G and 10"^^ G. 
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For the same two sets of parameters above these are given by the dashed and dotted curves 
respectively. The delayed spectra for 50 s duration (topmost dashed and dotted curves) 
corresponds to the GRB duration (tgrb) and we have calculated them using the full range of 
secondary pair energy 7e = e^/2me, created with CMB, allowed by the maximum primary 
photon energy ^ '-f'^^^^rrieTi) from equation (6). We have also calculated the delayed 
spectra for 10^, 10*^ and 10® s duration by numerically integrating equation (22) up to a 
maximum 7e corresponding to that time from Figure 3. 

5. Discussion 

We have calculated the spectrum and time dependence of the high energy (> GeV) 
spectrum emerging from GRB internal shocks as a function of the luminosity, time variability 
and the bulk Lorentz factor, using exact cross sections. For typical bursts parameters the 
internal shock spectrum cuts off above a threshold energy -E^,th ~ 10 — 100 GeV due to 
77 iteractions within the shock region. There is a qualitative difference with the results 
of Baring & Harding (1997), who discussed the 77 absorption in external shocks, and who 
did not consider the secondary spectra from interactions with background radiation. More 
recently Wang et al. (2004) also discussed the 77 absorption in external shocks, including 
interactions against background photons. The main difference with our work is due to the 
much lower 77 optical depth from the lower photon density at the radii of external shocks, 
which are much larger than the radii of internal shocks considered here. Internal shock 
77 absorption using approximate cross sections was discussed by Pilla & Loeb (1998) and 
Lithwick & Sari (2001) without interactions against external background photons, and by 
Dai & Lu (2002) including such interactions. These authors did not consider the effects of 
synchrotron self-absorption. The main qualitative difference between the emerging internal 
shock spectra of these authors and our work arises from our use of accurate cross sections, 
and more importantly, our inclusion of synchrotron self-absorption. The cut-offs in the 
emerging spectra that we obtain are compatible with more detailed numerical calculations 
of the self-consistent pair creation and annihilation spectrum in internal shocks of Pe'er 
& Waxman (2003), without inclusion of interactions against external background photons. 
The main difference between their final spectra and ours is due to their using an ab-initio 
spectrum instead of a phenomenological Band spectrum as here, and to the additional effects 
of our inclusion here of interactions againts external photons, leading to a delayed secondary 
spectrum. 

A new feature discussed in this paper is that the spectra that emerge from the GRB 
internal shock region shows, besides the expected high energy cutoff due to pair-production 
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£'-y,th, also a re-emergence of the spectrum at a higher energy E'^/^thin > -E'7,th where the shock 
becomes optically thin again to 77 interactions. This is due to our inclusion of synchrotron 
self-absorption, which reduces the internal photon target spectral density causing absorption 
of the highest energy emerging photons. For high enough bulk Lorentz factors (typically 
Tft ^ 800), the fireball is optically thin to internal 77 interactions at all energies, and there 
is no internal 77 cut-off band in the emerging spectrum. 

The emerging GRB spectrum is modified, on its way to the observer, by interactions with 
the diffuse infrared background (CIB) and the diffuse microwave background (CMB). Here, 
we have used the GRB internal shock emerging spectra including synchrotron self-absorption 
as the source of input photons, and calculate their interaction with a CIB spectrum rcdshifted 
to z ~ 1, based on that of Malkan & Stecker (2001). We then calculated the secondary GRB 
spectrum from the resulting secondary pairs which are up-scattered on CMB photons. We 
used two different escaping GRB primary spectra, one with an internal 77 cutoff band (for 
a typical burst bulk Lorentz factor Fb ~ 300) and one without an internal 77 cutoff band 
(for a burst with a larger Ff, ~ 800). We used two different isotropic-equivalcnt luminosities, 
and different intergalactic magnetic field strengths ranging between Biq = 10^^° G and 
10^^^ G. This leads to a secondary spectrum, mainly in the 1-100 GeV range, in addition to 
and delayed respect to the prompt unabsorbed internal shock spectrum. We also calculated 
the further interaction of the secondary gamma-ray spectrum with the CIB, which gives a 
similarly attenuated spectral shape as that of the primary spectrum. We did not consider 
the delayed spectrum of the tertiary and higher order pairs, since the delay time scales are 
longer and the fluxes of these components is much lower. 

A different delayed GeV component can arise in the afterglow phase, when the in- 
verse Compton spectrum peak energy in the external shock region sweeps across the GeV 
band (Meszaros & Rees 1994; Dermer, Chiang & Mitman 2000; Sari & Esin 2001; Zhang & 
Meszaros 2001). This delayed GeV emission, from IC upscattering on external shock elec- 
trons, is predicted to be detectable by GLAST at 2; ~ 1 (Zhang & Meszaros 2001), as long 
as the shock parameters are such that the IC component is prominent (which is the case 
for the typical parameters inferred from broadband modeling of some well-studied bursts 
(Wijers & Galama 1999; Panaitescu & Kumar 2002; Yost et al. 2003)). The IC duration 
can be up to a few hours, which could overlap with the earlier part of the delayed spectra 
discussed here. The temporal evolution of the IC spectral component is such that it hardens 
with time initially, but softens with time after the IC peak crosses the band. The delayed 
GeV component discussed in this paper has a rather different temporal evolution behavior. 
Since Af(7e) is essentially anti-correlated with 7e (see Figure 3), the hard spectrum (due 
to the more energetic electrons) lasts a shorter timescale than the soft spectrum. Also the 
harder photons arrive earlier than the softer photons. Throughout the whole delayed phase. 
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the hardest photons emerge first, followed by progressively softer photons. Later the hardest 
photons drop out first from the spectrum. The delayed spectrum always progressively softens 
with time. This distinct spectral evolution behavior can be used to distinguish the delayed 
emission from IC scattering in external shocks from the delayed emission of secondary spec- 
tra due to pair-production and IC scattering on the diffuse background discussed here. The 
secondary external spectrum from such an IC spectrum, however, should be similar to what 
we have discussed here, on timescales longer than a few hours, since the absorbed energy 
undergoes similar reprocessing. 

The time delays and the flux level of the secondary spectra arc sensitive to the inter- 
galactic magnetic field strength. The best prospects for detection of such secondary delayed 
GeV spectra occurs for bursts with bulk Lorentz factors > 800 and intergalactic magnetic 
fields < 10"^^ G, which might be typical of intergalactic void regions. Since the mean-free- 
path of the first generation high energy photons against 77 absorption typically exceeds 
several to tens of Mpc (Coppi & Aharonian 1997; Malkan & Stecker 2001), it is plausible 
that a substantial fraction of the pair formation will occur in void regions near the burst, 
even if the latter occured in a cluster galaxy. The corresponding time delays are in the 
range 10-100 s and the fluxes are within the sensitivity range of GLAST, and at the high 
energy end, also of VERITAS, HEGRA, and other ACTs. A detection of this secondary 
delayed GeV spectrum in GRB would provide a valuable diagnostic for the IGM magnetic 
field strength, as well as the bulk Lorentz factor, and would provide valuable clues for the 
typical size of the shocks responsible for the radiation. 

This research is supported in part by NSF AST 0098416, NASA NAG5-13286 and the 
Monell Foundation. 
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Fig. 1. — Opacity for high energy photons in the GRB fireball, in the co moving frame, 
against electron Compton scattering, pair production with electrons by 67 and by 77 
interactions with fireball photons. We have plotted the curves for L^^iso = 10^^ erg/s isotropic 
7-ray luminosity, the bulk Lorentz factor Tf, = 316 (solid curves), 500 (dotted curves) and 
800 (dashed curves). We used a variability time 6t = 0.01 s. Photons above ~ 2 x 10^ GeV, 
for example, for Tf, = 316 will escape the GRB fireball. No photons, for the same parameters, 
may escape the fireball below ~ 2 x 10'' GeV down to ~ 26 GeV where the 77 opacity is 
greater than unity. 
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Fig. 2. — Cartoon of allowed (optically thin) and forbidden (optically thick) photon energy 
range, in the observer's frame, from the GRB fireball as a function of the bulk Lorentz factor 
(Xb). Photons below -E^,th and above -E^^thin may escape. We have plotted the curves for 
2 different 7-ray luminosity 10^^L52, lO^^Lsi; 4 different variability time ISto s, 0.15t_i s, 
0.015t_2 s, 0.0015t_3 s; and 3 different peak synchrotron photon energy 20052.3 keV, 500e2.7 
keV, lOOO^s o keV. However, we have kept the synchrotron self-absorption energy to be fixed 
at 10 keV in all cases. The parameter choices are shown for curves from top to bottom. 
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Fig. 3. — Different time scales (Atjc, AIa, At^, and tgrb), involved in calculatng the delayed 
secondary IC spectrum using equation (22), as a function of secondary pair energy (7e). 
The thick gray lines correspond to At(7e) = max{Atic, At^, At^, St} for two different Ats, 
denoted by solid and dashed lines, corresponding to inter-galactic magnetic field strength 
(-Big) of 10~^^ and 10~^° G respectively. We have used a GRB duration of tgrb = 50 s. All 
times are plotted for a fixed luminosity distance of Di = 10^^ cm. 
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Fig. 4. — High energy photon spectrum {E^^^ = E'^(PN^/ dtdE^) from a GRB where the 
fireball is 77 optically thin to photons of all energy. The primary spectrum (full line) is due 
internal shocks, the secondary spectrum (dashed and dotted) is due to 77 interactions and 
IC scattering in the external environment. The parameters are 5t = 1 s, e^^pk = 1 MeV, = 
10^'^, isotropic-equivalent luminosity 10^^ erg/s, spectral index a = 2.5, located at redshift 
z = 1. The delayed spectra are calculated for two different intergalactic magnetic field values, 
BiQ = 10~^° and 10~^^ G, denoted by the dashed and the dotted curves respectively. The 
duration of the delayed secondary emission (from top to bottom) spectra are 50 s, 10^ s, 10^ 
s and 10^ s. The GLAST sensitivity (thick gray lines) for 50 s, 10^ s and 10^ s integration 
time is also plotted for comparison. 
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Fig. 5. — High energy photon spectrum {E^^^ = E'^cPN^/ dtdE^) from a GRB where the 
fireball is optically thick due to 77 against some its own photons over a certain energy band. 
The parameters are 5t = 0.1 s, e^^pk = 0.5 MeV, Ff, = 10^ '', isotropic-equivalent luminosity 
Xq52 erg/s, spectral index a = 2.5, and redshift z = 1. Plotted are the prompt (solid curve) 
and delayed spectra (dashed and dotted curves). The delayed spectra are for two different 
intergalactic magnetic fields of -Big = 10~^° and 10^^^ G, denoted by the dashed and the 
dotted curves respectively. The duration of delayed emission (from top to bottom) spectra 
are 50 s, 10^ s, 10^ s and 10^ s. GLAST sensitivities (thick gray lines) for 50 s, 10^ s and 
10^ s integration time are also plotted for comparison. 



